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基于 行星 际 闪 烁 (1IPS) 望远镜 的 数字 波束 
合成 技术 研究 
杨 轩 "， 陈 林 杰 '”， 荣 健 ， 磊 才 华 ” 
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北京 100049) 
摘 要 : 对 于 拟 建设 的 子午 工程 二 期 中 的 明 安 图 行星 际 闪烁 望远镜 ， 数 字 波 束 合成 技术 是 其 
中 的 一 项 关键 技术 。 针对 现 有 数字 波束 合成 算法 无 法 满足 IPS 望 远 镜 上 的 应 用 需求 和 复杂 度 较 高 ， 
本 文 对 现 有 的 数字 波束 合成 算法 进行 了 改进 和 优化 , 提出 了 一 种 适应 于 明 安 图 IPS 望 远 镜 的 波束 
合成 算法 ,仿真 结果 显示 该 算法 能 够 很 好 地 满足 IPS 望 远 镜 的 需求 ,同时 ,在 单个 波束 的 基础 上 ， 
本 文 也 提出 了 多 波束 合成 的 算法 ， 并 进行 了 算法 的 仿真 。 
关键 词 : 行星 际 闪烁 (IPS) 望远镜 ; ”数字 波束 合成 算法 ; ”线性 约束 最 小 方差 (LCMV) 算法 
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国家 重大 科研 装备 研制 项 目 子午 工程 二 期 列 入 了 国家 发 改 委 十 大 优先 支持 项 目 ， 该 工程 连 
续 监 测 地 球 表面 20km-30km 以 上 到 几 百 公里 的 高 层 大 气 、 电 离 层 和 磁 层 以 及 十 儿 个 地 球 半径 以 
外 的 行星 际 空间 环境 中 的 电磁 场 、 电 场 、 中 高 层 大 气 的 风 场 、 密 度 和 成 分 、 电 离 层 、 磁 层 和 行 
星际 空间 的 有 关 参 数 。 


行星 际 内 炮 (Interplanetary Scintillation: IPS) 是 指 来 自 于 遥远 致密 射电 源 的 电磁 波 在 穿 过 
OO 太阳 系 行星 际 空间 时 ， 被 太阳 风 等 离子 散射 形成 的 射电 波 强度 和 相位 的 一 种 随机 起 伏 现象 。 通 
pi 过 研究 这 种 闪烁 现象 ， 就 可 以 研究 致密 源 的 角 结 构 ， 并 测量 太阳 风速 度 和 太阳 风 等 离子 不 规则 
> 结构 ， 观 测 数 据 对 理论 研究 和 实际 应 用 都 有 重要 作用 *。 行 星际 闪烁 的 地 基 测 量 是 监测 行星 际 
空间 的 一 种 经 济 、 有 效 的 遥感 技术 。 以 天 然 射 电源 或 卫星 下 行 通 讯 链 为 射电 信 标 ， 以 地 基 望 远 
镜 的 单 站 或 多 站 组 网 方式 接收 射电 信号 ， 行 星际 空间 作为 射电 传播 路 径 的 介质 能 影响 射电 的 传 
播 特 性 。 通 过 对 该 透射 的 射电 的 信号 强度 、 传 播 时 延 、 相 位 变化 等 参量 的 测量 ， 必 要 时 加 以 反 
演 计算 ， 可 以 遥测 行星 际 空间 的 结构 及 其 时 空 变化 特征 。 
目前 ， 国 外 多 个 国家 已 经 利用 现 有 望远镜 开展 了 对 行星 际 闪 烁 的 研究 ， 国 内 由 中 科 院 空间 
中 心 牵头 建设 的 子午 二 期 行星 际 闪烁 观测 设备 目前 正在 进行 技术 储备 ， 其 通过 三 站 联合 观测 ， 
并 分 析 数 据 ， 有望 在 背景 太阳 风速 的 三 维 大 尺度 分 布 和 太阳 活动 周 演化 、 行 星际 CME 的 三 维 空间 
层 析 反 演 和 时 变 追 踪 、 太 阳 风 的 动 理学 油 流 幅度 等 方面 取得 重要 科学 进展 ”。 
对 于 行星 际 内 烁 GIPS) 望远镜 ， 其 中 的 一 个 关键 技术 是 数字 多 波束 合成 技术 。 


数字 多 波束 合成 亦 称 空域 滤波 器 ， 是 阵列 处 理 的 一 个 主要 技术 。 目 前 ， 它 在 射电 天 文 ， 雷 
达 、 声 呐 、 导 航 、 无 线 电 通信 和 地 震 监测 等 诸多 方面 都 得 到 了 广泛 的 应 用 。 本 文 主要 讨论 其 在 
射电 望远镜 上 的 应 用 。 数 字 多 波束 合成 实质 是 充分 利用 天 线 阵列 在 不 同位 置 的 阵 元 所 获得 的 空 
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间 信 号 信息 ， 通 过 对 各 阵 元 进行 加 权 ， 来 达到 增强 期 望 信号 、 抑 制 干扰 的 目的 。 其 优点 有 : 可 


术 可 形成 多 个 主 波 束 ， 实 现 对 天 空 多 个 射电 源 同时 进行 观测 。 
1 行星 际 闪 烁 CIPS) 望远镜 
在 子午 二 期 行星 际 闪烁 望 远 镜 的 设计 中 ， 分 析 比 较 日 本 Nagoya 四 站 


形成 单个 或 多 个 独立 可 控 的 波束 而 不 损失 信 品 比 ， 波 束 特性 由 权 矢 量 控制 ， 灵 活 多 变 ， 天线 具 
有 较 好 的 自 校正 和 低 旁 瓣 能 力 ;方便 后 续 进行 阵列 信号 处 理 ， 以 获得 优良 性 能 ”。 利 用 这 一 技 


和 印度 0oty 单 站 的 IPS 


设计 方案 和 观测 经 验 ， 我 们 提出 了 新 的 “一 主 站 两 辅 站 ”的 望远镜 设计 方案 。 其 主 站 位 于 内 蒙 


采用 全 可 动 的 抛物 面 天 线 。 该 三 站 的 地 理 坐 标 近 似 构 成 “等 边 三 角形 ”， 
长 度 近似 200 公 里 。IPS 探 测 设备 主要 是 测量 射电 源 强度 闪烁 的 时 间 序 列 ， 


古 正 镶 白 旗 明 安 图 镇 ， 采 用 大 型 可 转动 的 抛物 柱 面 天 线 ; 两 辅 站 分 别 位 于 阿 巴 嘎 和 克 什 殉 腾 ， 


其 任 两 站 之 间 的 基线 
采用 单 台 站 观测 的 功 


率 谱 拟 合 和 三 台 站 联 测 的 互相 关 时 延 ， 快 速 获 得 不 同 空 间 区 域 处 太阳 风 的 速度 、 电 子 密度 扰动 
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的 非 约 化 功率 谱 等 物理 参数 ， 追 踪 太 阳 风 和 CME 扰 动 在 行星 际 空 间 的 加 速 


图 1 IPS 望 远 镜 的 布局 
Fig.1 Configuration of the IPS Telescope 


、 传 播 和 演变 。 


子午 二 期 ITPS 望 远 镜 三 个 站 点 中 的 主 站 由 三 个 抛物 柱 面 天 线 并 排放 置 
(N-S) 方向 上 的 尺寸 为 140 米 ， 在 东西 〈E-W) 方向 上 的 尺寸 为 40 米 ， 如 


图 2 IPS 望 远 镜 主 站 阵列 
Fig.2 main site array of IPS Telescope 


扫描 天 空 和 寻找 其 他 科学 目标 ， 柱 面 天 线 被 设计 成 可 以 在 E-W 方 向 转动 ， 


, 每 个 柱 面 天 线 在 南北 
图 2 所 示 。 为 了 快速 


但 不 跟踪 天 空中 的 源 。 


通过 转动 ， 在 不 同 的 时 间 指 向 不 同 的 天 区 ， 可 以 完成 整个 天 区 的 快速 扫描 覆盖 。 另 外 ， 在 主 站 
还 有 一 个 20 米 的 抛物 面 天 线 ， 也 可 以 用 于 IPS 的 观测 。 两 个 辅 站 各 有 一 个 约 20 米 的 抛物 面 天 线 。 


在 每 个 站 点 ， 天 线 不 仅 可 以 双 频 率 单 站 点 独立 地 观测 射电 源 ， 而 且 也 可 以 双 频 率 多 站 点 联合 工 
作 。 通 过 这 种 方式 ， 一 方面 可 以 利用 主 站 的 天 线 获得 很 高 的 灵敏 度 ， 观 测 到 天 裤 中 更 多 的 射电 
源 ， 另 - 方面 还 可 以 利用 对 强 射电 源 进 行 多 站 相关 来 获得 更 为 精确 的 太阳 风速 度 等 信息 ， 从 而 
更 好 地 重建 三 维 太阳 风 结 构 。 
每 个 柱 面 天 线 有 一 个 i 沿 着 案 焦 线 放 置 304 个 馈 源 接收 射电 信号 ， 可 以 在 327MHz 和 
654MHz 下 同时 或 是 单独 工作 。 每 个 馈 源 连接 着 一 个 低 噪声 放大 器 LNA ) ， 信 号 经 放大 、 带 通 
滤波 后 经 过 一 个 AD 转换 器 进行 数字 化 。 数 字 信和 号 通过 光纤 传输 到 室内 接收 机 ， 然 后 进行 数字 多 
波束 合成 ， 最 后 输出 不 同方 向 多 个 射电 源 的 观测 信和 号。 本 文 主要 是 对 IPS 望 远 镜 主 站 柱 面 天 线 的 


数字 多 波束 合成 技术 进行 研究 ”。IPS 望 远 镜 阵列 的 主要 参数 和 技术 指标 如 表 1 所 示 。 
表 1 1PS 望 远 镜 阵列 的 主要 参数 和 技术 指标 


Tab. 1 Main Specifications and technical indicators of 1PS Telescope array 


天 线 柱 面 天 线 : 3X140mX40m 
抛物 面 天 线 : 1X20m 和 2X“20m 

频率 327MHz, 654MHz, 1420MHz” 

频率 带宽 40MHZ 

频率 分 辨 率 39KHZ 

合成 波束 主 站 : 3*12 个 ， 旁 办 低 于 -20dB 

灵敏 度 主 站 :“4. 4mJy ， 子 站 : “261mJy 


(9) 仅 辅 站 抛物 面 天 线 在 此 频率 观测 
2 阵列 信号 的 建 模 

假设 空间 有 一 个 射电 源 ， 且 有 N 个 阵 元 组 成 的 均匀 线 阵 阵列 来 接收 信号 。 将 第 一 个 阵 元 位 
于 坐标 原点 ,作为 参考 阵 元 。 在 信号 是 窄带 的 情况 下 , 若 参 考 阵 元 接收 的 信号 用 SG) = s(t)eY 表 
示 ， 则 第 i 个 阵 元 接收 的 信号 可 以 表示 为 5,(1) = s(t 一 t,)ewe W'。 均 匀 线 阵 接收 信号 模型 如 区 


| 


ti=dsin®@ /c 


CO 


Eee 


2 


4 所 示 。 


图 4 均匀 线 阵 接收 信号 模型 


Fig.4 Received signal model of uniform linear array 


i j —jwt — jwt 二 
我 们 设 a(O) 为 导向 矢量 ， 则 4a(O) =| Le ee se | ,其 中 第 个 阵 
元 相对 于 参考 阵 元 的 延 时 上 =id sin 0/c 。 


干扰 和 期 望 信号 为 同 频率 的 罕 带 信号 ， 但 入 射 角度 不 同 ， 且 与 期 望 信号 不 相关 。 


噪声 信号 为 平稳 的 、 零 均值 ,方差 为 o” 的 高 斯 自 噪声 , 并 且 各 阵 元 接收 的 噪声 间 互 不 相关 ， 
且 于 期 望 信号 不 相关 。 


3 阵列 方向 图 定义 


阵列 天 线 的 方向 图 : 是 按 一 定 方式 排列 的 天 线 阵 列 的 权 矢 量 对 来 自 空间 不 同方 向 的 信号 的 
阵列 响应 。 则 阵列 的 方向 函数 可 以 表示 为 
F(0)= W”a(O) (1) 
式 中 W 为 波束 合成 算法 得 出 的 权 矢 量 。 
方向 函数 表征 了 波束 合成 后 ， 天 线 对 各 个 方向 来 的 信号 的 增益 。 对 上 式 进 行 取 绝 对 
值 并 归 一 化 ， 然 后 取 对 数 处 理 后 ， 阵 列 的 方向 图 可 以 表示 为 


P(O=20log 二 (2) 


max|F (0)| 
由 () 和 (2) 式 可 以 得 出 ， 当 阵列 的 权 矢 量 改 变 时 ， 对 来 自 不 同方 向 的 空间 信和 号 来 说 , 就 会 
有 不 同 的 响应 ， 从 而 形成 的 阵列 方向 图 P(O) 的 指向 就 会 改变 。 


4 数字 波束 合成 算法 


呆 数字 波束 合成 算法 是 通过 对 阵列 不 同位 置 的 阵 元 进行 加 权 ， 来 使 阵列 方向 图 的 主办 对 准 期 
望 信号 入 射 方向 ， 低 旁 闪 对 准 干扰 方向 ， 以 此 达到 增强 期 望 信 和 号、 抑制 干扰 的 目的 ， 其 包括 党 
后 规 数 字 波束 合成 算法 以 及 自 适应 数字 波束 合成 算法 。 


> 4.1 常规 数字 波束 合成 算法 

j 常规 数字 波束 合成 算法 的 权 值 是 固定 的 , 与 输入 信号 无 关 。 权 值 等 于 期 望 信 号 的 导向 矢量 。 
因此 在 期 望 信号 方向 的 响应 为 F(6@,) = a(0@,)”a(0,) = ;其 中 ,为 期 望 信号 入 射 角 , N 为 阵 元 
个 数 。 这 时 ， 各 通道 的 信号 同 相 相 加 ， 阵 列 在 期 望 信号 0 方向 的 信号 增益 最 大 。 

4. 2 自 适应 数字 波束 合成 算法 

自 适 应 数字 波束 合成 算法 是 根据 信号 的 先 验 知识 和 一 定 准 则 ， 调 节 阵列 的 权 矢 量 ， 使 天 线 
= 阵列 方向 图 在 期 望 方向 形成 主 波束 ， 而 在 干扰 方向 形成 零 陷 ， 从 而 对 干扰 形成 抑制 。 自 适应 波 
© 束 合成 权 值 是 与 输入 信号 有 关 的 ， 且 随 着 电磁 环境 的 变化 而 自 适应 的 调整 ， 能 够 实现 对 期 望 信 
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号 的 跟踪 “"”。 自 适应 波束 合成 系统 框图 如 图 5 所 示 。 
图 5 自 适应 波束 合成 系统 框图 


Fig.5 Block diagram of adaptive beamforming system 
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LCMYV 算法 是 一 种 较 优 的 自 适应 数字 波束 合成 算法 , LCMYV 算法 原理 是 在 控制 期 望 信号 方 
向 的 增益 为 1， 干扰 方向 的 增益 为 0 的 前 提 下 ， 使 阵列 输出 功率 最 小 外 。 其 约束 表达 式 为 


C”W= 1 (3) 


式 中 C=[a(@),a(Q)] 为 Nx2 维 约束 和 矩阵， 为 期 望 信 号 方向 ，0 为 干扰 信号 方向 ， 
f=[1,0] 为 2x1 维 响 应 矩阵 。 
LCMYV 算法 的 数学 表达 为 


min(W” RW) (4) 


SC 了 =f (5) 


式 中 民 , 为 接收 信号 (期 望 信号 、 干 扰 和 噪声 ) 的 协 方差 矩阵 。 
利用 拉 格 衣 日 乘 子 法 解 得 最 优 权 矢 量 关 


W,,=R CIC"R AC| / (6) 


由 于 系统 幅 相 误差 的 存在 ， 会 使 接收 到 的 信号 的 协 方差 矩阵 尺 ,和 实际 信号 的 协 方差 矩阵 
R, 存在 差异 ， 算 法 的 性 能 也 会 因此 下 降 。R, 的 误差 也 就 是 期 望 信号 的 误差 和 干扰 的 误差 ， 为 
了 减少 误差 的 影响 , 对 R, 做 一 个 预 处 理 , 我 们 知道 LCMV 算法 是 在 期 望 信号 响应 为 1 的 情况 下 ， 
让 期 望 信号 、 干扰、 噪声 的 功率 最 小 ， 也 就 是 让 干扰 和 噪声 的 功率 最 小 。 所 以 ， 在 RR, 中 减 去 
望 信号 的 成 分 


X =X.-X. (7) 


R=E(X, X2 ) (8) 


tn itn 


式 中 义 , 为 接收 到 的 信号 (期望 信 号 、 干 扰 和 噪声 ) ，X. 为 接收 到 的 期 望 信号 ，X,,, 为 接收 到 的 
eg 


Wo=R, CIC"R,C| f (9) 


opt0 ~ Titn 


这 样 接收 到 的 信号 的 协 方差 矩阵 RR 和 实际 信号 的 协 方差 矩阵 R, 的 误差 就 是 接收 到 的 干扰 
和 实际 的 干扰 之 间 的 误差 。 
4. 3 Matlab 实 验 仿真 与 结果 分 析 

对 于 IPS 望 远 镜 主 站 的 抛物 柱 面 天 线 , 其 单个 柱 面 天 线 包含 304 个 馈 源 , 馈 源 间 距 为 半 波 长 。 
在 IPS 观 测 中 ， 需 要 对 所 有 馈 源 的 信号 进行 波束 合成 处 理 ， 来 实现 天 线 指向 不 同 的 射电 源 。 这 里 
可 以 将 所 有 这 些 馈 源 视 为 一 个 馈 源 线 阵 ， 阵 元 个 数 N=304 个 ， 阵 元 间距 d=0.5 入 。 为 了 评估 数字 
波束 合成 算法 的 性 能 ， 我 们 进行 了 如 下 仿真 : 
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(1) 假定 期 望 信号 到 达 角 为 10 度 。 常 规 数字 波束 合成 算法 下 的 阵列 方向 图 如 图 6 所 示 。 
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图 6 常规 波束 合成 算法 下 形成 的 一 个 主 波束 


Fig.6 A main beam formed by conventional beamforming algorithms 


(2) 期 望 信号 到 达 角 为 10 度 ， 干 扰 信 号 到 达 角 为 50 度 ， 快 拍 数 为 1000，SNR=30dB, INR=30dB， 


Amplitude in dB 


-20 


Amplitude in dB 
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干扰 加 上 5% 的 幅 相 误差 。 基 于 LCMV 算 法 下 的 阵列 方向 
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图 如 图 7 所 示 。 


图 7 LCMV 算 法 下 形成 的 一 个 主 波束 
Fig.7 A main beam formed by LCMV _ algorithms 

仿真 结果 显示 LCMV 算 法 可 以 抑制 干扰 ， 而 常规 波束 合成 算法 不 能 抑制 干扰 。 两 种 算法 阵 
列 方向 图 的 第 一 旁 六 都 较 高 ， 达 到 -13.4dB 左 右 ， 不 满足 IPS 望 远 镜 数字 波束 合成 旁 洲 -20dB 的 技 
术 指 标 要 求 〈 表 1) 。 在 观测 中 ， 干 扰 有 可 能 从 较 高 的 旁 流 进入 接收 系统 ， 有 影 响 观测 结果 。 
5 改进 的 数字 多 波束 合成 算法 

常规 波束 合成 算法 可 以 通过 加 窗 的 方法 降低 旁 办 以 达到 IPS 望 远 镜 技 术 指 标 ， 但 加 窗 的 方 
法 损失 了 天 线 的 灵敏 度 。 本 文 在 LCMV 算 法 的 基础 上 提出 了 一 种 改进 算法 ， 在 不 损失 天 线 灵 敏 
度 的 情况 下 降低 旁 汶 ， 以 满足 IPS 望 远 镜 技 术 指 标的 要 求 。 
5. 1 改进 的 LCMV 算 法 

由 LCMV 算 法 可 知 ， 其 可 以 在 干扰 信号 方向 形成 深 零 陷 。 为 此 ， 采 用 设置 虚拟 干扰 方法 使 
劳 流 降 低 。 在 LCMV 算 法 结果 上 ， 和 寻找 高 于 某 个 值 (比如 -25dB) 旁 久 峰 值 所 在 的 角度 ， 并 将 其 设 
为 虚拟 干扰 信号 方向 LCMV 算法 是 让 干扰 信号 的 响应 为 0， 即 虚拟 干扰 方向 深 零 陷 ， 但 这 必 
使 其 他 方向 位 置 有 较 大 升 高 。 所 以 ， 改 进 算法 让 虚拟 干扰 响应 不 为 0， 而 为 一 个 很 小 的 数 〈 比 
H0.05) ， 这 样 即 可 以 让 虚拟 干扰 方向 的 旁 流 高 度 降低 ， 也 不 会 让 虚拟 干扰 附近 的 部 分 旁 流 高 
有 较 大 升 高 。 那 么 改进 的 LCMV 算 法 的 约束 表达 式 可 以 表示 为 

C”"W=/f (10) 

式 中 C =[a(@,),a(0),.….,a(0,)] 为 Nx(m+l) 维 约束 矩阵 ，@ 位 期 望 信号 方向 ，@ 为 干扰 信号 
方向 ,多 …@, 为 虚拟 干扰 方向 。 卢 = 岂 ,0,0.05,……,0.05] 为 (m+1)xl 维 响应 矩阵 。 

改进 的 LCMYV 算法 的 数学 表示 为 


六 


阔 呈 


Ek 


通 


minC W” RW) (11) 
stC"W=f (12) 

利用 拉 格 朗 日 乘 子 法 解 得 最 优 权 矢量 为 
WR OR | 站 (13) 


5. 2 Matlab 实验 仿真 与 结果 分 析 

在 同样 的 条 件 下 ， 假 定期 望 信号 到 达 角 为 10 度 ， 干 扰 信 号 到 达 角 为 50 度 。 快 拍 数 为 1000， 
SNR=30dB, INR=30dB, 干扰 加 上 5% 的 幅 相 误差 。 将 高 于 -25db 的 角度 设置 为 虚拟 干扰 方向 。 采 
用 改进 的 LCMV 算 法 ,仿真 结果 如 图 8 所 示 。 
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图 8 改进 的 LCMV 算 法 下 形成 的 一 个 主 波束 
Fig.8 A main beam formed by improved LCMV algorithms 


20 


仿真 采用 了 改进 的 LCMV 算 法 ,将 LCMV 算 法 下 的 阵列 方向 图 中 高 于 -25dB 的 角度 设置 为 虚 
拟 干 扰 方向 ， 并 让 其 响应 为 0.05。 仿 真 结 果 可 以 看 到 阵列 方向 图 主办 附近 的 高 于 -25dB 的 劳 准 高 
度 降低 至 -253dB 以 下 。 由 于 增加 了 虚拟 干扰 方向 响应 的 约束 条 件 ， 虚 拟 干 扰 方 向 附近 的 部 分 旁 兴 
的 幅度 有 小 幅度 的 增加 ， 使 其 超过 -25dB， 但 没有 超过 -20 dB。 因 此 ， 改 进 后 的 LCMV 算 法 达到 


了 IPS 望 远 镜 的 技术 指标 (-20dB) 要 求 。 
6 基于 子 阵 的 改进 数字 波束 合成 算法 


如 果 对 阵列 的 每 个 阵 元 都 对 应 一 个 接收 通道 ， 这 时 阵列 的 响应 最 好 ， 阵 列 控制 波束 指向 的 
灵活 度 也 最 高 。 但 是 ， 对 于 本 文 这 种 大 型 阵列 权 矢 量 的 维度 就 会 很 高 ， 进 而 数字 波束 合成 算法 
的 复杂 度 就 会 很 高 ， 同 时 每 个 接收 通道 都 需要 连接 AD 转换 器 、 滤 波 器 等 便 件 资源 ， 系 统 非常 复 
杂 ， 便 件 成 本 也 很 高 。 于 是 希望 尽 可 能 保持 良好 的 阵列 响应 前 提 下 ， 通 过 减少 通道 数 的 方式 来 
降低 系统 复杂 度 和 硬件 成 本 。 减 少 通道 数 是 采用 子 阵 化 方法 实现 的 ， 每 个 子 阵 对 应 一 个 接收 通 


道 ”。 因 此 ， 下 面 提出 基于 子 阵 的 数字 波束 合成 算法 。 基 于 子 阵 的 数字 多 波束 合成 算法 原 到 
图 如 图 9 所 示 。 


E 框 


阵 元 
移 相 器 
子 阵 
| 
+ 


求 和 输出 


| 


图 9 基于 子 阵 的 数字 多 波束 合成 算法 原理 框图 


Fig.9 Principle Block diagram of digital multi-beam beamforming algorithm based on subarray 


表 2 阵 元 级 和 子 阵 级 复杂 度 对 比 


原始 阵列 子 阵 方案 | (4 个 阵 元 ) 子 阵 方案 11 (8 个 阵 元 ) 
阵列 单元 数 304 76 38 
模拟 移 相 器 0 76X3 38X7 
AD 转换 器 304 76 38 
数据 速率 304 X 100MB 76 X 100MB 38X 100MB 


权 矢 量 维度 304X1 76X1 38X1 


Tab. 2 Comparisons of complexity at element level and subarray Level 


上 表 是 阵 元 级 和 子 阵 级 复杂 度 对 比 ， 可 以 看 到 采用 基于 子 阵 的 数字 波束 合成 算法 ， 除 了 子 
阵 内 阵 元 之 间 需 要 加 多 个 模拟 移 相 器 外 ， 在 AD 转 换 器 个 数 、 数 据 速 率 、 权 矢量 维度 都 以 单个 子 
阵 内 阵 元 个 数 成 倍 的 降低 。 单 个 子 阵 内 的 阵 元 数目 越 多 ，AD 转 换 器 个 数 、 数 据 速 率 、 权 矢量 维 
度 降 低 的 越 多 ， 但 是 阵列 控制 波束 的 灵活 度 又 有 所 降低 。 所 以 为 了 保持 展 好 的 阵列 响应 ， 不 能 
一 味 的 增加 子 阵 内 阵 元 的 个 数 。 


6. 1 Matlab 实 验 仿真 与 结果 分 析 


(1) 对 于 由 304 个 阵 元 构成 的 线 阵 ， 假 定子 阵 阵 元 数 分 别 为 19、38， 子 阵 间 无 重 登 ， 期 望 信号 到 
达 角 为 10 度 ， 干 扰 到 达 角 为 50 度 。 快 拍 数 为 1000，SNR=30dB, INR=30dB， 干 扰 加 上 5% 的 幅 相 
误差 。 将 高 于 -25dB 的 角度 设置 为 虚拟 干扰 方向 。 阵 元 之 间 采 用 移 相 器 ， 使 同一 子 阵 内 阵 元 相位 
为 同 相 。 子 阵 之 间 选 用 改进 的 LCMV 算 法 。 对 子 阵 内 不 同 阵 元 数 分 别 进行 仿真 ， 仿 真 结果 如 图 
10、11 所 示 。 
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图 10 基于 子 阵 (19 阵 元 的 改进 LCMV 算 法 下 形成 的 一 个 主 波束 
Fig. 10 A main beam formed by improved LCMV algorithm based on subarray (19 elements) 
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图 11 基于 子 阵 (38 阵 元 ) 的 改进 LCMV 算法 下 形成 的 一 个 主 波束 
Fig.11 A main beam formed by improved LCMV algorithm based on subarray (38 elements) 
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仿真 结果 显示 基于 子 阵 (单个 子 阵 内 19 阵 元 和 38 阵 元 ) 的 改进 LCMV 算 法 都 满足 了 IPS 望 远 镜 
的 旁 欠 指标 要 求 ， 算 法 的 复杂 度 相 比 阵 元 级 数字 波束 合成 算法 也 大 大 降低 。 但 是 ， 以 38 个 阵 元 
为 一 个 子 阵 的 波束 合成 算法 在 干扰 方向 并 没有 形成 深 霉 陷 ， 且 第 一 旁 欠 的 幅度 比较 高 ， 算 法 的 
性 能 明显 降低 。 还 对 基于 子 阵 ( 子 阵 阵 元 数 分 别 为 2、4、8、16) 的 改进 LCMV 算 法 进行 了 仿真 ， 仿 
真 结果 也 都 满足 了 IPS 望 远 镜 的 劳 闪 指 标 要 求 。 而 且 ， 随 着 单个 子 阵 阵 元 数 的 增加 ， 算 法 性 能 会 
越 来 越 低 。 因 此 在 可 以 整除 304 的 数字 中 ，19 为 最 大 的 单个 子 阵 阵 元 数 。 


7 多 波束 合成 

前 面 主要 讨论 了 单个 波束 的 合成 ， 而 我 们 需要 同时 对 多 个 射电 源 进 行 观测 ， 就 需要 形成 多 
个 主 波束 。 要 形成 多 个 主 波束 ， 根 据 我 们 前 面 所 述 的 波束 合成 原理 可 知 ， 阵 列 的 权 值 决定 了 阵 
列 方向 图 主 波束 的 波束 指向 。 因 此 ， 就 需要 产生 多 组 不 同 的 权 值 。 每 组 权 值 对 应 不 同 的 相位 补 


Wn Wm Wn 
求 和 输出 | 求 和 输出 
ly C(t) | ym C(t) 


偿 和 幅度 调整 值 ， 从 而 形成 多 个 主 波 束 " 。 合 成 多 波束 的 原理 框图 如 图 12 所 示 。 
图 12 多 波束 合成 的 原理 框图 


Fig.12 Principle block diagram of multi-beam beamforming 


7. 1 Matlab 实 验 仿真 与 结果 分 析 


(1) 仿真 条 件 : 阵 元 个 数 N=304 个 ， 阵 元 间距 d=0. 5 入 ，19 个 阵 元 为 一 个 子 阵 无 重 着 ， 多 个 期 望 
言 号 到 达 角 为 分 别 为 -20 度 、0 度 、10 度 、30 度 , 干扰 信号 到 达 角 为 50 度 。 快 拍 数 为 1000, SNR=30dB， 
INR=30dB， 干扰 加 上 5% 的 幅 相 误差 。 将 高 于 -25db 的 角度 设置 为 虚拟 干扰 方向 。 阵 元 之 间 采 用 移 
相 器 ， 使 同一 子 阵 内 阵 元 相位 为 同 相 ， 子 阵 之 间 选 用 改进 的 LCMV 算 法 。 仿 真 结果 如 图 13 所 示 。 
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图 13 基于 子 阵 〈19 阵 元 ) 的 改进 LCMV 算 法 下 形成 的 4 个 主 波束 
Fig.13 Four main beams formed by improved LCMV algorithm based on sub-array (19 elements) 

由 上 图 可 知 ， 在 -20 度 、0 度 、10 度 、30 度 4 个 方向 形成 了 主 波束 ， 且 劳 兴 也 基本 达到 了 IPS 
望远镜 技术 指标 。 
8 总 结 

本 文 对 子午 二 期 IPS 望 远 镜 的 关键 技术 一 数字 波束 合成 技术 进行 了 研究 , 提出 了 一 种 适应 于 
IPS 望 远 镜 的 柱 面 天 线 的 数字 波束 合成 算法 。 该 算法 在 LCMV 算 法 的 基础 上 利用 设置 虚拟 干扰 的 方 
法 获得 了 满足 了 IPS 望 远 镜 的 合成 波束 。 同 时 ， 采 用 子 阵 化 的 方法 降低 了 系统 数字 波束 合成 算法 
的 复杂 度 , 进而 降低 了 工程 实现 难度 。 最 后 , 在 单个 波束 的 基础 上 , 提出 了 多 波束 合成 的 算法 ， 
实现 对 天 空 多 个 射电 源 同 时 进行 观测 。 
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Abstract: For the Mingantu interplanetary Scintillation telescope of the Meridian II project the 
digital beamforming technology is one of the key technologies. Since the present digital beamforming 
algorithms cannot match the requirements of the interplanetary Scintillation telescope, In this paper, 
the existing digital beamforming algorithm - LCMYV algorithm is Improved and optimized. The 
method of setting virtual interference is used to reduce the sidelobe of array pattern and the method of 
subarray is used to reduce the complexity of the algorithm. The improved LCMYV algorithm is a digital 
beamforming algorithm suitable for Mingantu interplanetary scintillation telescope. The simulation 
results show that the algorithm can meet the needs of IPS telescope very well. Meanwhile, based on 
the single beam algorithm, an algorithm of multi-beamforming has also been proposed and discussed. 


Key words: Interplanetary scintillation (IPS) telescope, Digital beamforming algorithm, 
Linear Constrained Minimum Variance (LCMYV) 


